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We determine the curvature of the pseudo-critical line of strong interactions by means of numerical 
simulations at imaginary chemical potentials. We consider Nf = 2 -|- 1 stout improved staggered 
fermions with physical quark masses and the tree level Symanzik gauge action, and explore four 
different sets of lattice spacings, corresponding to Nt = 6,8,10,12, in order to extrapolate results 
to the continuum limit. Our final estimate is k = 0.0135(20). 
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I. INTRODUCTION 


The exploration of the phase diagram of strongly in¬ 
teracting matter in the temperature - baryon chemical 
potential (T —/is) plane is being pursued both by exper¬ 
imental and by theoretical investigations. The compar¬ 
ison between the chemical freeze-out line and the 

crossover line, corresponding to chiral symmetry restora¬ 
tion, is one of the main issues. In principle these two lines 
are not expected to coincide, however an exact statement 
about their interrelation will provide useful information 
about the dynamics of strong interactions. That requires 
a precise determination of both lines. 

From the theoretical point of view, Lattice QCD sim¬ 
ulations represent the best first principle tool to provide 
information about the chiral transition^ temperature T^: 
present results provide consistent evidence for Tc cs 155 
MeV at pLB = 0. Unfortunately, as one moves to finite 
baryon chemical potential, direct numerical simulations 
are presently hindered by the so-called sign problem, 
stemming from the complex nature of the fermion de¬ 
terminant when fiB ^ 0. However, various methods have 
been proposed to circumvent the problem in the regime 
of small chemical potentials, where the pseudo-critical 
line can be well approximated by a quadratic behavior^ 
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^ We speak of chiral transition even if present lattice studies pro¬ 
vide evidence for a crossover [^13||. 

^ We note that a possible ambiguity in the denominator of the 
quadratic term, i.e. whether we take /Tc{^b) oi" /^s/ 77(0) as 
an expansion variable, is irrelevant as long as just the quadratic 
term is considered, since it only affects higher order terms. 


in 

Tc{fb) ^ 4 , 

—— = ( 1 ) 

where the coefficient k defines the curvature of the 
pseudo-critical line Tc{fJ-B)- Information about n can 
be obtained for instance by Taylor expansion tech¬ 
niques EMI): by analytic continuation from imagi¬ 
nary chemical potentials [isl - l^ . by reweighting tech¬ 
niques [sol or W a reconstruction of the canonical 
partition function [sj, . 

Recent numerical investigations [2^, , adopting the 

method of analytic continuation with improved dis¬ 
cretizations at or close to the physical point of Nf = 2-1-1 
QCD, have provided results for n which are generally 
larger than previous estimates obtained by the Taylor 
expansion technique [15 41711 . 

In particular, in Ref. 29| we performed numerical sim¬ 
ulations adopting an improved stout staggered fermion 
discretization on lattices with Nt = 6, 8, leading to a pre¬ 
liminary estimate n ^ 0.013, to be compared with previ- 
ous determinations obtained by Taylor expansion |l5l4l7j| , 
reporting k ^ 0.006. 

In the present study we aim at extending our results 
in two directions. First, we increase the number of 
imaginary chemical potentials explored on lattices with 
Nt = 8, in order to obtain a better control over the an¬ 
alytic continuation systematics and to perform a deeper 
comparison between the cases in which a strange quark 
chemical potential is included or not. Then we extend 
simulations for /is = 0 to two new sets of lattice spac¬ 
ings, corresponding to Nt = 10 and Nt = 12, in order 
to perform a continuum extrapolation of our determina¬ 
tion of K. As a byproduct, we also discuss the behavior 
of the continuum extrapolated chiral susceptibilities as a 
function of /is, in order to assess the possible influence 
of the baryon chemical potential on the strength of the 
transition, which is relevant to the possible existence of 
a critical endpoint in the T — pi b plane. 

The paper is organized as follows. In Section [IT] we pro¬ 
vide some details about the lattice discretization adopted 
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in this study, about the various explored setups of chem¬ 
ical potentials, about the observables chosen to locate 
and their renormalization. In Section IlIIl we discuss our 
numerical results and finally, in Section lIVi we draw our 
conclusions. 


II. NUMERICAL SETUP 


As in Ref. [2^, we consider a lattice discretization of 
Nf = 2-1-1 QCD in the presence of purely imaginary 
quark chemical potentials. We consider the following eu¬ 
clidean partition function 


Z = 


f—u, d, s 




Sym 

(Mi),,, = 


0. - IJ' 1J I- 


,3 + 0 


(4) 


where U are the gauge link variables, Sym is the tree 
level improved Symanzik gauge action [s^ , written 
in terms of (trace of the n x m loop constructed 

from the gauge links along the directions v depart¬ 
ing from the i site). Finally, the staggered Dirac opera¬ 
tor is built up in terms of the two times stout- 

smeared [ 3 ^ links with an isotropic smearing pa¬ 

rameter p = 0.15. Stout smearing improvement is used in 
order to reduce taste symmetry violations (see Ref. [s^] 
for a comparison among different improved staggered dis¬ 
cretizations) ; the rooting procedure is exploited, as usual, 
to remove the residual fourth degenerany of the staggered 
lattice Dirac operator (see, e.g.. Ref. Q for a discussion 
on possible related systematics). 

The temperature of the system is given by T = 
l/{Nta), where a is the lattice spacings and Nt is 
the number of lattice sites in the temporal direction, 
along which we take thermal boundary conditions (pe- 
riodic/antiperiodic for boson/fermion fields). At fixed 
Nt, T is changed by varying the value of the bare cou¬ 
pling constant /3. The bare quark masses mg and mi are 
rescaled accordingly, in order to move on a line of con¬ 
stant physics, with ttItt — 135 MeV and mg/mi = 28.15. 
This line is determined by a spline interpolation of the 
values reported in Refs. [H, 113 (see also Ref. [^1. 
Four different sets of lattice spacings, corresponding to 
Nt = 6, 8,10,12, have been explored, in order to extrap¬ 
olate our results to the continuum limit. 


A. Setup of chemical potentials 

In Eq. ([2]), we have introduced an imaginary chemical 
potential pf = ipjj, Pfj G M, with / = u,d,s, cou¬ 
pled to the number operator of each quark flavor. They 
are related to the chemical potentials coupled to con¬ 
served charges (baryon number B, electric charge Q and 
strangeness S) by the following relations 

l+u = pb/^ + ^pq/S 

= pB/i-^J,Q/?> ( 5 ) 

Ps = pB/i - Mq/3 - PS ■ 

The purpose of our study is to determine the dependence 
of the pseudocritical temperature Tc on the baryon chem¬ 
ical potential (which is given hy pB = Pu + ‘2pd), in a 
setup of chemical potentials which is as close as possible 
to the thermal equilibrium conditions created in heavy 
ion collisions. We thus have to require to S' = 0 and 
Q = rB, where r is the number of protons divided by the 
number of nucleons of the colliding ions, r = ZjA k, 0.4 
typically. 

These requirements can be translated into relations 
between pB, ps and pQ, which at the lowest order in 
Pb read pQ ~ qi(T)pB and ps — Si{T)pB, the coef¬ 
ficients qi{T) and si(r) being related to derivatives of 
the free energy density [4l|, [dg. Let us consider as an 
example the strangeness neutrality condition: in a gas 
of non-interacting fermions it would imply ps = 0 but 
in QCD, due to interactions, the mixed derivatives of 
the free energy density with respect to ps and pu, Pd are 
non-vanishing, so that one needs a non-zero p,,, to en¬ 
sure 5 = 0. Present lattice investigations 1^ show 
that, for T ^ 155 MeV, the constraints on charge and 
strangeness imply si ~ 0.25 and qi ~ —0.025. With a 
precision of a few percent, around the transition at van¬ 
ishing density, we thus have pi = pu = Pd, Pi — Pb/3 
and Ps ^ Pi/4:. 

Our determination of the curvature k has been ob¬ 
tained setting Ps = 0, which is close to the conditions 
described above. To quantify the impact of ps, as in 
Ref. [13 , we have considered also the case ps = pi, in or¬ 
der to obtain an estimate about the effect of a non-zero 
Ps in a range which covers the equilibrium conditions 
created in heavy ion collisions. 


B. Physical observables, renormalization and the 
determination of Tc 

In the absence of a true phase transition, the deter¬ 
mination of the pseudo-critical line may depend on the 
physical observable chosen to locate it. On the other 
hand, chiral symmetry restoration is the leading phe¬ 
nomenon around Tc, with the light chiral condensate be¬ 
coming an exact order parameter in limit of zero light 
quark masses. Therefore in the following Tc{pb) will 
be determined by monitoring the chiral properties of the 
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system. The chiral condensate of the flavor / is defined 
as 


W)/ 


TdXogZ 
V druf 


( 6 ) 


where V is the spatial volume. In our simulations the 
two light quarks are degenerate, m/ = and it 

is convenient to introduce the light quark condensate: 


is affected by both additive and multiplicative 
renormalizations. We consider two different renormal¬ 
ization prescriptions, in order to determine whether any 
systematic effect related to this choice affects the deter¬ 
mination of K. The first one [1^ is 


W)(i)(T) 



- w(S'S) 

(T) 


(T = 0) 


( 8 ) 


where ms is the bare strange quark mass; in this way 
the leading mass dependent contribution is subtracted^, 
while one takes care of the multiplicative renormaliza¬ 
tion by dividing by the same quantity at T = 0. As an 
alternative, we consider the following prescription [l^ 


operator corresponding to a single light flavor) 




disc 


dA.SC _j_ com 

dmi ~ 


( 10 ) 


T fNi 
V 


[((TrMf 1)2) - (TrMf 1)2] (11) 


L-C// C I L 


T Ni , 

--^(TrM-). 


( 12 ) 


where Ni = 2 is the number of degenerate light quarks. 
The renormalization is performed by first subtracting the 
T = 0 contribution, to remove the additive renormaliza¬ 
tion, then multiplying the result by the square of the bare 
light quark mass, to cancel the multiplicative one [s^: 


[x^^T) - x^^{T = 0)] . (13) 


C. Analytic continuation from imaginary chemical 
potentials 

The physical observables relevant to our study will be 
monitored as a function of T for fixed values of the di¬ 
mensionless ratio 9i = Im(/r;)/T. In this way we shall 
be able to locate Tc for a set of values of 9i, so as to 
determine the dependence Tc{9i) to the leading order 

^^ = l + Ref+0{9f), (14) 


(00)(2) = ^ ((0V')z - {A'^)i(T = 0)) . (9) 

In this case the zero T subtraction eliminates additive 
divergences while multiplication by the bare quark mass 
mi takes care of multiplicative ones. 

The behavior of both condensates will be monitored to 
locate Tc- In particular, since in the presence of a true 
phase transition the slope of the condensate as a function 
of T diverges at T^, we will look for the point of maximum 
slope, i.e. the inflection point (a detailed comparison with 
other prescriptions has been reported in Ref. Hi). 

A much better probe is provided by the chiral suscepti¬ 
bility Xipipj which is itself divergent at in the presence 
of a true transition: in this case the introduction of rel¬ 
evant parameters (finite mass or finite volume) smooths 
the divergence, however looking for the maximum of Xip^ 
remains a well defined and univoque prescription for lo¬ 
cating the pseudo-critical temperature Tc- On the lattice, 
the light chiral susceptibility is given by (M; is the Dirac 


® This prescription subtracts both divergent and finite terms which 
are linear in the mass, thus permitting to isolate contributions 
to the quark condensate directly related to spontaneous chiral 
symmetry breaking. However, possibile additive logarithmic di¬ 
vergences could still be present. 


where we have assumed Tc{9i) to be an analytic function 
of 9i, at least for small values of it. This assumption 
is consistent with numerical data and is at the basis of 
the method of analytic continuation. Comparing with 
Eq. ID) one has, at the leading order in k = R/9. 

III. NUMERICAL RESULTS 

We have performed simulations on lattices with Nt = 
8 ,10 and 12 and different choices of T and of the chem¬ 
ical potentials; results will be combined with those al¬ 
ready presented in Ref. [1^ for Nt = 6,8 to perform the 
continuum extrapolation. To that purpose, we will con¬ 
sider only lattices with fixed aspect ratio Lg/Nt = 4: that 
guarantees the absence of significant finite size effects (see 
Ref. for a detailed study about that). 

Four different values of chemical potentials have been 
considered for Nt = 10,12, corresponding to /r® = 0 and 
Im(^/)/(7rT) = 0,0.20,0.24 and 0.275. A larger set has 
been considered for Nt = 8, in which case we performed 
simulations also at A 0, in order to provide more 
information about systematics related to the choice of 
[isl[i-i and to the truncation of the Taylor expansion in 
Eq. (flil) . 

For each setup of chemical potentials we have explored 
0(10) different temperatures around Tj9i). The Ratio¬ 
nal Hybrid Monte-Carlo algorithm |4hl - l4^ has been used 
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FIG. 1: Renormalized susceptibility and chiral condensates for the 40^ x 10 (left column) and 48® x 12 lattices (right colnmn). 


to sample gauge configurations according to Eq. each 
single run consisting of 2-5 K trajectories of unit length in 
molecular dynamics time, with higher statistics around 
the transition. 

Traces appearing in the definition of chiral quantities 
(see, e.g., Eqs. (ITTI) and (IT^ ) have been computed by 
noisy estimators at the end of each molecular dynamics 
trajectory, using 8 random vectors for each flavor. Such 
a choice has appeared, after some preliminary tests, as 
a reasonable compromise to balance the effort spent in 
the stochastic estimators and in the gauge configuration 
production, i.e. in order to optimize the statistical error 
obtained for a given computational effort. A jackknife 


analysis has been exploited to determine the statistical 
errors. 

To perform the renormalization described in Sec. m 
one needs to compute observables also at T = 0 and 
at the same values of the bare parameters, i.e. at the 
same ultraviolet (UV) cutoff. Eor that reason we have 
performed simulations on lattices as large as 48^: details 
are reported in Appendix 

In order to determine the inflection point of the renor¬ 
malized chiral condensate, we have performed a best fit 
to our data according to 

= Ai + Bi arctan(C'i(r - TJ) , (15) 











5 


Lattice 

Mi,/ 

ttT 

ttT 

Tc{'ipip(i)) 


Tcixl 

16® 

X 

6 

0.00 

0.00 

148.2(3) 

148.4(4) 

150.7(4) 

16® 

X 

6 

0.20 

0.00 

155.0(4) 

155.1(5) 

157.0(4) 

16® 

X 

6 

0.24 

0.00 

158.9(4) 

159.1(4) 

160.0(4) 

16® 

X 

6 

0.275 

0.00 

161.2(4) 

161.5(4) 

162.7(4) 


X 

6 

0.00 

0.00 

149.0(6) 

149.0(6) 

151.6(5) 

24® 

X 

6 

0.24 

0.00 

160.8(7) 

160.7(5) 

162.0(5) 

24® 

X 

6 

0.275 

0.00 

164.1(4) 

164.3(5) 

165.9(4) 


X 

6 

0.00 

0.00 

149.1(7) 

149.4(4) 

152.0(4) 

32® 

X 

6 

0.24 

0.00 

160.2(3) 

160.4(2) 

162.7(4) 

32® 

X 

6 

0.275 

0.00 

163.4(3) 

163.5(3) 

165.5(4) 

“sF 

X 

8 

0.00 

0.00 

154.2(4) 

154.5(4) 

155.6(7) 

32® 

X 

8 

0.10 

0.00 

155.4(7) 

155.2(8) 

157.2(7) 

32® 

X 

8 

0.15 

0.00 

159.5(9) 

158.9(9) 

160.2(5) 

32® 

X 

8 

0.20 

0.00 

162.9(8) 

163.0(6) 

163.0(6) 

32® 

X 

8 

0.24 

0.00 

165.0(5) 

164.8(5) 

165.8(8) 

32® 

X 

8 

0.275 

0.00 

169.5(9) 

168.6(7) 

169.8(7) 

32® 

X 

8 

0.30 

0.00 

172.4(9) 

171.8(9) 

172.8(8) 

32® 

X 

8 

0.10 

0.10 

157.1(8) 

157.0(8) 

158.5(7) 

32® 

X 

8 

0.15 

0.15 

159.2(9) 

158.8(8) 

160.1(8) 

32® 

X 

8 

0.20 

0.20 

163.9(6) 

163.7(6) 

165.3(9) 

32® 

X 

8 

0.24 

0.24 

169.4(7) 

168.6(6) 

169.6(7) 

32® 

X 

8 

0.275 

0.275 

175.4(6) 

174.4(7) 

177.0(8) 

40®” 

X 

10 

0.00 

0.00 

154.5(1.5) 

154.3(1.5) 

155.1(7) 

40® 

X 

10 

0.20 

0.00 

163.0(7) 

163.0(8) 

162.5(7) 

40® 

X 

10 

0.24 

0.00 

166.8(8) 

167.1(7) 

166.2(1.0) 

40® 

X 

10 

0.275 

0.00 

170.8(8) 

171.2(8) 

169.6(8) 

48®” 

X 

12 

0.00 

0.00 

154.5(1.0) 

155.5(1.3) 

154.7(7) 

48® 

X 

12 

0.20 

0.00 

163.2(1.2) 

165.0(1.5) 

161.9(7) 

48® 

X 

12 

0.24 

0.00 

165.2(1.1) 

166.2(1.0) 

166.2(1.0) 

48® 

X 

12 

0.275 

0.00 

167.8(1.2) 

168.7(9) 

167.9(9) 


TABLE I: Critical values of T obtained from the renormalized 
chiral susceptibility and from the renormalized chiral conden¬ 
sates. Errors do not take into account the uncertainty on the 
physical scale, which is of the order of 2-3% [3^. |40|. 

which involves the independent parameters Ai, Bi, Ci 
and Tc- Instead, for the determination of the peak of the 
renormalized susceptibility, we have performed a best fit 
according to a Lorentzian function 

= 52 (T - Tcf ■ 

Both functions are found to well describe respective data 
points around T^. In both cases, statistical errors on 
the fitted parameters have been estimated by means 
of a bootstrap analysis, while systematic uncertainties 
have been estimated either by varying the range of fit¬ 
ted points around Tc or by choosing an alternative fitting 
function (e.g., a hyperbolic tangent for the condensate or 
a parabola for its susceptibility). Statistical and system¬ 
atic^ errors are both included in the collection of deter¬ 
minations of Tc for the various combinations of lattice 


We do not report the systematic error on the determination of the 



1 ' 1 ' 1 

0 

^ 1 ^ 1 

- n <viiv|t>p, 


0 X 


I 


j: 







_ 

I 

• 

_1_,_1_,_1 

_ ^^^^_ 


J_^^^_I_^_I_^_L_ 

0 0.02 0.04 0.06 0.08 

(p,/(3tT))" 


FIG. 2: Critical lines for the 40^ x 10 lattice (top) and for the 
48® X 12 one (bottom). 


sizes and chemical potentials in Table HI which includes, 
for completeness, also results presented in Ref. [^ . 

In Fig. [T] we report results obtained for ('i/''0)(i) 
and ('0V')(2) the 40® x 10 and 48® x 12 lattice, together 
with some best fits according to Eqs. m and (ITSl) . In 
the following we will perform the continuum limit using 
two different methods, in order to check for systematics 
effects. 


Lattice 

«('!/’ V’(l)) 

K(t/)V’(2)) 

«(x’') 

24® X 6 
32® X 8 
40® X 10 
48® X 12 

0.0150(7) 

0.0142(7) 

0.0157(17) 

0.0130(15) 

0.00152(7) 

0.0135(7) 

0.0164(16) 

0.0123(17) 

0.0140(7) 

0.0134(9) 

0.0139(10) 

0.0131(11) 


TABLE IL Curvatures obtained at fixed Nt from different 
observables. 


physical scale, which is of the order of 2-3 % [ 3 ^ . Iioll and, being 
related to an overall scale determination, does not affect the ratio 
of pseudocritical temperatures entering the determination of k, 
see Eqs. 0 and m. 
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FIG. 3: Continuum limit of the curvatures extracted at fixed 
Nt (data have been slightly shifted in the horizontal direction 
to improve readability). 

A. Continuum limit for ps = 0 - First method 

In order to extract the curvature of the critical line, 
we have performed a best fit to the values of 
obtained for each lattice size and setup of chemical po¬ 
tentials, according to the function 

Ufiij) = TM + ■ 

(17) 

For all sets of chemical potentials explored for = 0, the 
inclusion of quartic corrections has not been necessary: 
a more detailed discussion about the stability of the fit 
as the range of chemical potentials is changed is reported 
in Sec. IIIIDI 

In Fig. [5] we report an example of such quadratic fits 
to the critical temperatures obtained for Nt = 10,12 and 
for the various explored observables. A complete collec¬ 
tion of results, including also those already presented in 
Ref. [13, is reported in Table HIl 

In a range of temperatures around Tc, the UV cut¬ 
off a~^ is approximately proportional to Nt- Therefore, 
assuming corrections proportional to a^, we extracted, 
from the curvatures obtained for different values of Nt, 
continuum extrapolated results according to the ansatz 

K{Nt) = Kcont + const./Nt . (18) 

Results are shown in Fig. [31 where we also re¬ 
port the extrapolated continuum values, which 
are «:cont((^V')[i)) = 0.0134(13), Kco™t((V'^/^)( 2 )) = 
0.0127(14) and K,out{x/j,^) = 0.0132(10). 

B. Continuum limit for ps = 0 - Second method 

Results of the previous section show that the contin¬ 
uum extrapolation of k is quite smooth, with a good 


agreement between the results obtained with different 
observables and different renormalization prescriptions. 
This is also consistent with the preliminary evidence re¬ 
ported in Ref. [^ . 

Nevertheless, it is useful to explore different ways of 
performing the continuum limit, in order to check for 
the overall consistency of the procedure. In the previous 
section we first determined the value of k at each single 
value of Nt, then extrapolated these results to ^ oo to 
obtain Kcont- A different procedure is to first extrapolate 
the critical temperatures to iVt —>■ oo (for fixed values 
of the dimensionless ratio ^ij/T) and then to extract 
the value of Kcont by using the continuum extrapolated 
critical temperatures. 

To implement the second procedure we have per¬ 
formed, separately for each ^ij/T, a best ht to the val¬ 
ues obtained for the renormalized condensates and for 
the renormalized chiral susceptibility on different values 
of Nt, according to modified versions of Eqs. (ESI and 
EH- Since the cut-off dependence is more pronounced 
for such quantities, we have excluded Nt = 6 data, thus 
using only Nt = 8,10,12. 

In detail, in the case of the renormalized susceptibil¬ 
ity, each fit parameter appearing in Eq. EH has been 
given an additional Nt dependence, for instance Tc{Nt) = 
Tc{Nt = oo) -I- const/N(^. Results for the extrapolated 
quantities are reported in the upper plot in Eig.|3] where, 
for the sake of clarity, we report only the cases fiij = 0 
and fj,ij/{TrT) = 0.275. In the case of the renormalized 
condensates, instead, due to the larger number of param¬ 
eters which are present in Eq. ESI. we could obtain fits 
which are stable against the variation of the fitted range 
by adding an A^(-dependence to just two parameters, in 
particular Tc and Ci. Results are shown in the middle 
and lower plot of Eig. 01 

Such fits provide estimates for the continuum extrapo¬ 
lated pseudo-critical temperatures, reported in Table m 
and in Eig. |SJ Such values coincide, within errors, 
with the continuum pseudo-critical temperatures that 
one could obtain by directly fitting results reported in Ta- 
bleH A best fit to the extrapolated temperatures accord¬ 
ing to Eq. EZI, with only the quadratic term included, 
provides Kcont((^'0)(1)) = 0.0145(11), Kcont(('0V')(2)) = 
0.0138(10) and Kcont{x/p^) = 0.0131(12), which are con¬ 
sistent with those found previously. 


fJ.l.l/i'KT) 

To(V>V’(i)) 

Tc(V'V'(2)) 

W) 

0.00 

0.20 

0.24 

0.275 

154.7(8) 

163.9(8) 

166.9(9) 

169.7(8) 

156.5(8) 

165.0(7) 

168.5(7) 

170.8(7) 

154.4(8) 

161.0(1.1) 

165.8(1.0) 

167.3(1.1) 


TABLE III: Continuum extrapolated critical temperatures for 
the various pi j values. 
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0.02 - 










Com. exi. |i^^/(7rT) = 0.0 
Com. ext. p.^^/(7rT) = 0.275 

0 Nt=12 
o Nt=10 
□ Nt = 8 
Nt= 12 
Nl= 10 
Nt = 8 




140 160 


, , 

180 200 
T [MeV] 


0.8 - 


0.4 - 


Com. ext. |j.^^/(7iT) = 0.0 



160 180 
T [MeV] 



160 180 
T [MeV] 


180 

175 

170 

3 165 
160 

155 


0 0.02 0.04 0.06 0.08 

(H„/(31T))" 

FIG. 5: Critical lines obtained by using the continuum extrap¬ 
olated renormalized chiral susceptibility and the continuum 
extrapolated chiral condensates. 



pect a visible dependence of the strength parameters also 
for small values of imaginary ^b- The width and the 
height would tend respectively to zero and infinity ap¬ 
proaching, e.g., a critical endpoint in the Z 2 universality 
class. 

To that purpose, in Fig. [6] we plot the continuum ex¬ 
trapolated width B 2 and height A 2 /BI as a function of 
fiij. No apparent change of either quantity can be appre¬ 
ciated, hence no dependence of the strength as a function 
of /rs. 

Of course, that does not exclude the presence of a crit¬ 
ical endpoint at real fiB- the critical region could be 
small enough, or the endpoint location far enough from 
^B = 0, so that no influence is visible for small, imagi¬ 
nary ^B- For instance, for = 0, a Roberge-Weiss 
like endpoint is expected along the pseudo-critical line at 
imaginary chemical potential, for ^i j/{nT) ^ 0.45 [ 2 ^. 
Fig. [5] shows that also this endpoint has no apparent in¬ 
fluence on the strength of the transition in the explored 
range. 


FIG. 4: Gontinuum limit for the renormalized susceptibility Yy Inclusion of /i, / 0 and systematics of analytic 
and the renormalized chiral condensates. continuation 


C. Strength of the transition as a function of fiB 

The width and the height of the chiral susceptibility 
peak, which can be obtained respectively from B 2 and 
A 2 /B 2 in Eq. (fTBl) . are directly related to the strength 
of the chiral pseudo-transition. Therefore, we have the 
possibility to monitor the dependence of such strength on 
the baryon chemical potential and, having performed a 
continuum extrapolation for can do that directly 

on continuum extrapolated quantities. 

If a critical endpoint exists, along the pseudo-critical 
line, for relatively small values of real /i^, we might ex- 


We have extended results for Nt = 8 presented in 
Ref. [^ . performing numerical simulations for a larger 
range of imaginary chemical potentials, which include 
also the case /is = /i;. That enable us to answer two 
important questions. What is the systematic error, in 
the determination of k by analytic continuation, related 
to the truncation of the Taylor series in Eq. (II3 and to 
the chosen range of chemical potentials? What is the im¬ 
pact of our effective ignorance about the actual value of 
/is corresponding to the thermal equilibrium conditions? 
We are going to discuss in detail only the determination 
of the pseudo-critical temperature from the renormal¬ 
ized chiral susceptibility, however we stress that similar 
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FIG. 6: Peak values (xlOO) and widths of the continuum 
extrapolated renormalized chiral susceptibility. 

conclusions are reached when one considers the renor¬ 
malized chiral condensate. The corresponding pseudo- 
critical temperatures, taken from Table HI are reported 
in Fig. [7] for fig = 0 and for = it/. 

We first tried a quadratic fit in remembering the 
defintion 0/ = nij/T, we used 

Tc(0/) =r,(O)(l + 9K0f) (19) 

and several fits have been performed by changing each 
time the maximum value included in the fit. Rea¬ 

sonable best fits are obtained in all cases, apart from the 
fit to the whole its = M/ range, which yields a reduced 
^ 2.4 and indicates the need for quartic corrections 
in this case. Results obtained for k are shown in Fig. |S1 
for Us = 0, the fitted value of k is perfectly stable as 
the range of chemical potentials is changed. Instead, for 
Its = It/, the value of tt clearly depends on the fitted range 
of chemical potentials: it is larger as the range is ex¬ 
tended and becomes compatible, within errors, with that 
obtained for its = 0 as the range is decreased. This be¬ 
havior is consistent with the presence of significant quar¬ 
tic corrections in this case. That may be related to the 
different structures of the phase diagrams for imaginary 
chemical potential that one has in the two cases: this 
issue has been discussed in detail in Ref. [^ . 

We then tried a best fit to a function including quartic 
corrections, 

T,{0i) = T,iO)il + 9K0f+ bet), ( 20 ) 

to the whole range of chemical potentials explored in both 
cases. The corresponding results obtained for k are re¬ 
ported in Fig. m as well. While for 1 /^ = 0 the value 
is perfectly compatible with the one obtained without 
including quartic corrections (indeed, in this case one 
obtains b = 0 within errors), for itg = it/ we observe 
a significant change, bringing k in good agreement with 
the Its = 0 case. A similar conclusion is reached when 


^ o 32^x8 g, = 0 
° 32^x8 = 


J 





FIG. 7: Gritical lines for the 32® x8 lattices in the two different 
setups: Its = 0 and fia = fii- 

a common fit to both sets of data (i.e. with a common 
value for 113 ( 0 )) is performed, as shown in the right panel 
of Fig. IS] and in Fig. [T] 

We conclude that, for its = 0, no evidence of quar¬ 
tic corrections is found in the whole explored range. As 
a consequence, the extracted k is stable against varia¬ 
tions of the fitted range and we can exclude the presence 
of significant systematic corrections, related to the pro¬ 
cedure of analytic continuation, affecting the continuum 
extrapolated determination of k that we have provided. 

In the case its = li/, larger values of k are obtained 
when quartic corrections are neglected, however k be¬ 
comes compatible with that obtained for its = 0 when 
such corrections are included, or when the fitted range of 
chemical potentials is small enough. We conclude that 
K is not affected by the inclusion of its, at least within 
present errors, which however are larger than for the 
Its = 0 case. In particular, a fair estimate in this case is 
k{hs = M/) = 0.013(3). 


IV. CONCLUSIONS 

In the present study, we have extended results reported 
in Ref. [29| by performing numerical simulations on lat¬ 
tices with Nt = 10,12 and aspect ratio 4, and by enlarg¬ 
ing the range of chemical potentials explored for Nt = 8. 
That has permitted us to obtain continuum extrapolated 
results and to better estimate possible systematics re¬ 
lated to analytic continuation. 

Regarding the case fis = 0, we have obtained contin¬ 
uum extrapolated values of k from different observables 
(chiral susceptibility and the chiral condensate with two 
different renormalization prescriptions) and by two dif¬ 
ferent extrapolation procedures (extrapolating Kcont from 
k(W) or extracting Kcont from continuum extrapolated 
temperatures). The comparison of the two different pro¬ 
cedures permits us to give an estimate of the system¬ 
atic uncertainties related to the continuum extrapola- 
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0.04403(5) 

3.55 

32^ 

1.97(5) 

0.05680(13) 

0.03164(7) 

3.60 

32^^ 

2.05(6) 

0.03912(14) 

0.02211(7) 

3.65 

32^ 

1.82(3) 

0.02633(2) 

0.01518(9) 

3.70 

32^^ 

1.80(3) 

0.01804(3) 

0.01064(2) 

3.65 

48'‘ 

1.74(7) 

0.02638(4) 

0.01521(2) 

3.75 

48^^ 

1.61(5) 

0.01232(5) 

0.00749(2) 

3.85 

48"^ 

1.47(4) 

0.00614(2) 

0.00401(1) 

3.95 

48^^ 

1.37(3) 

0.00331(2) 

0.00237(1) 


TABLE IV: Determination of the observables at T = 0 (on the 
32^ and 48"^ lattices) needed to perform the renormalizations 
discussed in Section im Data are in lattice units. 


FIG. 8: Stability analysis of the ht to extract the k value 
for the 32 X 8 lattice. Empty symbols correspond to purely 
quadratic fit while filled symbols also take into account the 
quartic correction; red circles represents the = /r; data, 
black squares the Hs = 0 ones. The right panel shows the 
result of a combined ht (i.e. hxing a common value for Tc{0)) 
to both data sets when a quartic correction is used for the 
fia = fJ^i data: the open (hlled) triangle corresponds to fia = 0 

(/Ts = m). 


tion. In the case of the renormalized chiral susceptibility 
(k = 0.0132(10) vs K = 0.0131(12)) the systematic er¬ 
ror is negligible in comparison to the statistical one. In 
the case of ('0'i/')(i) = 0.0134(13) vs k = 0.0145(11)) 

and of (k = 0.0127(14) vs k = 0.0138(10)) 

the systematic and statistical uncertainties are clearly 
comparable in size. The extended analysis performed 
on iVt = 8 has permitted us to state also that, within 
present errors, systematic effects connected to the range 
of /x; chosen to extract the curvature are not signifi¬ 
cant. Regarding finite size effects, the analysis reported 
in Ref. already showed that they are negligible within 
the present precision on lattices with aspect ratio 4. Tak¬ 
ing into account the obtained results and the contribu¬ 
tions from the systematic effects mentioned above, we 
quote K = 0.0135(15) as our hnal continuum estimate for 
the case = 0. 

Such a result confirms, even after continuum extrap¬ 
olation, a discrepancy with previous determinations ob¬ 
tained by Taylor expansion [isl - ll^ . reportiM k ~ 0.006. 
As already discussed quantitatively in Ref. [29| , only part 
of this discrepancy can be accounted for by the different 
prescriptions used to determine the dependence of on 
/i;. Contrary to the Taylor expansion case, when work¬ 
ing at imaginary /r; one can use consistently the same 
prescription to locate Tc used for /i; = 0, i.e. looking 
for the maximum of the chiral susceptibility or the in¬ 
flection point of the chiral condensate (see Ref. [2^ for 
more details). The remaining part of the discrepancy 
could be possibly attributed to the systematic uncertain¬ 
ties related to the continuum extrapolation of previous 
studies. However, we stress that updated investigations 


by the same groups lead to results which are consistent 
with our estimate (see, e.g., Ref. [48l|). 

Regarding the case = /i;, we have confirmed the pre¬ 
liminary results reported in Ref. [2^. There is evidence 
for the presence of quartic contributions in the depen¬ 
dence of Tc on the imaginary in this case and when 
such contributions are taken into account, or when the 
range of htted chemical potentials around fj-B = 0 is small 
enough, the curvature becomes compatible, even if within 
larger errors, with that obtained for = 0. That means 
that also for the equilibrium conditions created in heavy 
ion collisions, corresponding to /ig 0.25 /i; around Tc, 
one does not expect significant deviations from the results 
obtained for /Xg = 0: a prudential estimate for the cur¬ 
vature in this case is^ k = 0.0135(20). That is obtained 
based on the estimate for /Xg = 0, with an increased error 
determined on the basis of the uncertainty that we have 
for the curvature extracted at /Xg = /x;. 

Finally, the analysis of the continuum extrapolated 
peak of the chiral susceptibility as a function of imag¬ 
inary ^B shows no significant varations of the strength 
of the transition, which could be associated to a possible 
nearby critical endpoint present along the pseudo-critical 
line. 
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Appendix A: Data at T = 0 

The determination of the renormalized condensate 
and susceptibility requires the computation of the cor¬ 
responding quantities at T = 0 and at the same UV 
cutoff of the finite temperature data. To that aim, we 
spanned a range of /3 on the line of constant physics, 
3.5 < /? < 3.95. The lattice sizes have been chosen in 
such a way to have temperatures well below Tc, keeping 
at the same time finite size effects under control. This 
required us to perform simulations on larger lattices (go¬ 
ing from 32^ up to 48^) as we decreased the value of the 
lattice spacing. We report results in Table HVl 
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FIG. 9: Comparison between the T = 0 and T ^ 0 chiral 
susceptibility ~ 0. The T = 0 susceptibility is needed 
to compute the renormalized chiral susceptibility, Eq. Ca}. 
Data are in lattice units and a linear fit to the T = 0 data is 
shown. 


The temperatures, which are in the range T ^ 25 — 50 
MeV, are low enough to be considered as a good approx¬ 
imation of the T = 0 limit; indeed, as expected because 
of the absence of transitions in this T range, observables 
depend smoothly on /3; moreover no dependence at all 
is expected on the imaginary chemical potentials, since 
they can be viewed as a modihcation in the temporal 
boundary conditions which, at T = 0 (i.e. for infinite 
temporal extension), are completely irrelevant. Hence, 
the relatively coarse sampling of the interval is enough 
to permit a reliable interpolation. We adopted a cubic 
spline interpolation for the condensate and a linar fit for 
the susceptibility. 


The renormalization prescription for the susceptibility, 
Eq. (fT^ . requires the subtraction of the T = 0 result from 
the hnite T contribution. To give an idea of the relative 
magnitude of this subtraction, in Fig. [9] we plot Xip-ip for 
zero chemical potential and both at zero and finite T. 
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